Alternating Kinetics of Annihilating Random Walks Near a Free Interface 
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The kinetics of annihilating random walks in one dimension, with the half-line x > initially 
filled, is investigated. The survival probability of the n th particle from the interface exhibits power- 
law decay, S n (i) ~ t~ a " , with a n ~ 0.225 for n = 1 and all odd values of n; for all n even, a faster 
decay with a n w 0.865 is observed. From consideration of the eventual survival probability in a 
finite cluster of particles, the rigorous bound ai < 1/4 is derived, while a heuristic argument gives 
Qi « 3\/3/8 = 0.2067. . .. Numerically, this latter value appears to be a lower bound for ai. The 
average position of the first particle moves to the right approximately as 1.7 t 1 / 2 , with a relatively 
sharp and asymmetric probability distribution. 



PACS numbers: 02.50.Ga, 05.70.Ln, 05.40. +j 

Annihilating random walks (ARWs) represent a sim- 
ple but ubiquitous reaction process in which particles dif- 
fuse and annihilate whenever they meet In addition 
to providing general insights about non-equilibrium phe- 
nomena, ARWs underlie a variety of basic kinetic pro- 
cesses ranging from the voter model M and the kinetic 
Ising-Glauber model u] , to reaction-diffusion systems Q 
and wetting phenomena In one dimension, power- 
ful exact solution methods have been developed to un- 
derstand many kinetic and spatial properties of ARWs 

While much is known about ARWs under homogeneous 
conditions, the role of spatial heterogeneity in such non- 
equilibrium systems is relatively unexplored. Our par- 
ticular interest is to understand the influence of a free 
interface on the asymptotic behavior of ARWs. For equi- 
librium systems at criticality, the presence of such a free 
interface gives rise to well understood surface critical be- 
havior which is characterized by associated surface criti- 
cal exponents |]l2|,^3). For reactive systems, in contrast, 
while there has been some progress in understanding the 
role of heterogeneity in the intrinsic properties of the re- 
actants , the role of heterogeneity in the form of 

a free interface is still unexplored. 

In this paper, we investigate basic properties of a one- 
dimensional semi- infinite system of ARWs. We consider a 
linear chain in which one particle initially occupies each 
lattice site n for n > 0, while the system is empty for 
n < 0. Far from the interface, the behavior should co- 
incide with that of the homogeneous system. However, 
the interfacial region is a less reactive environment be- 
cause one side of the system is initially empty. One 
might thus anticipate that particles near the interface 
should exhibit slower kinetics and different spatial prop- 
erties than bulk particles. This expectation is only partly 
correct; in fact, every other particle near the interface 
exhibits faster kinetics compared to bulk particles. Our 
work further indicates that there are only two apparently 
independent "surface" exponents which characterize the 
asymptotic particle survival probabilities. This surface 
behavior eventually governs the entire system, although 



it penetrates slowly into the bulk by diffusion. 

To quantify the phenomena that are governed by the 
existence of the interface in the semi-infinite ARW sys- 
tem, our work is organized around the following basic 
questions: 

• What is the probability that the first particle sur- 
vives until time t, Si(f)? More generally, what is 

,th 



the survival probability for the 
the interface, S n (t)? 



n 



particle from 



• What is the probability that particles i and j re- 
act as a function of \i — j\ (with i = 1 or 2 and j 
arbitrary)? 

• What is the spatial density distribution near the 
interface? 

To answer the first question, let us introduce the ex- 
ponents a n to characterize the probability that the n th 
particle survives until time t, S n (t) ~ t~° n . We first 
argue that there are only two independent exponents - 
one for n odd and a second for n even. To support this 
assertion, it is instructive to examine finite particle sys- 
tems. For 3 particles, the first and third particles have 
a finite probability to survive indefinitely, while the sec- 
ond certainly dies, with a survival probability which de- 
cays as t~ 3 / 2 [ ^0|p9|j20[ | . Similar asymptotic behavior can 
be anticipated for all finite systems with an odd num- 
ber of particles. In such cases, the asymptotically dom- 
inant contribution to S n (t) for n even will come from 
3-particle configurations where an even particle is be- 
tween two odd particles, with all other particles already 
annihilated. Conversely, particles with odd labels have 
a finite probability to survive indefinitely. Thus for the 
finite-particle system a n = for n odd and a n = 3/2 for 
n even. On this basis, we anticipate that just two ex- 
ponents also characterize the individual particle survival 
probabilities in the semi-infinite system. 

To test these predictions, we performed numerical sim- 
ulations using two complementary methods. The first 
(naive) approach is to simulate a suitably-sized system 
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with the initial condition that the right half-line is com- 
pletely occupied while the left half-line is empty. Re- 
flecting boundary conditions are employed at the edges 
of the system. The system size is chosen such that the 
effect of the boundaries is negligible over the time scale 
of the simulation. In the second approach, particles are 
created at the right boundary at a rate equal to the ex- 
act time-dependent density of the homogeneous system, 
c(t) ~ (47rt) -1 / 2 [ppc|], to mimic the effect of a semi- 
infinite system. This is a more efficient approach, as 
relatively long-time simulations can be run on a small 
systems without being influenced by boundary effects. 
For the survival probability of the first particle our sim- 
ulations give Si(t) ~ t~ ai , with a x = 0.225 ± 0.005. On 
the other hand, for the second particle, Sz(t) ~ t~ a2 with 
CK2 = 0.865 ± 0.015. As anticipated, the survival prob- 
ability of the first particle (n = 1) decays more slowly 
than t -1 / 2 , the particle survival probability in the bulk. 
However, the second particle is much less likely to survive 
than a bulk particle. This arises because the second par- 
ticle always has a unique potential left reaction partner, 
as well as a right reaction partner. 
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FIG. 1. Simulation results for Si(t) (o) and S 2 {t) (A). The 
quoted exponents are based on the best-fit straight lines to 
the data given in the figure. 

To determine whether the survival probabilities S n (t) 
are characterized by only two exponents, we perform a 
scaling analysis. Particles far from the interface initially 
exhibit bulk behavior, where S n (t) ~ t~ x l 2 for all n. Af- 
ter a time t n ~ n 2 , the n th particle "senses" the interface, 
and crossover from bulk to surface kinetics should occur. 
Based on the observed asymptotic behavior of Si(t) and 
S2(t), together with the above crossover picture, we ex- 
pect that S n (t) should exhibit the two distinct scaling 
forms for odd and even n respectively, 

S 2 „-i(t) * t" 1/2 0(«), S 2n (t) ~ t~ 1 ' 2 £ (z), (1) 

for t — > oo and n — > oo, with z = nt~ x l 2 finite. Large 
values of z corresponds to particles sufficiently deep in 



the bulk that they are not yet influenced by the het- 
erogeneous initial condition. Thus the large argument 
behavior of the scaling functions is determined by the 
survival probability of the homogeneous system. Since 
S(t) = c(t) ~ (4 7 ri)- 1 /2 |Q j this gives ^(oo) = 

£ (00) = 1/V47T- Conversely, for z -C 1 

00) £(z)~z^ 7 (2) 

with \i\ = 2a\ — 1 and /i2 = 2c*2 — 1 to match with the 
long time behavior of S\ (t) and S2 (t) . As Figure 2 shows, 
S n (t) follows this general description. 







A 




A 




D 




■ \ 

'is 




% 


An 

°° *%fa>ii?P*ligi]\ [ iiBi a 






m 





' 0.0 1.0 2.0 3.0 4.0 5.0 



FIG. 2. Simulation results for the time dependence of the 
survival probability of the n th particle, S n {t), for 1 < n < 100. 
Shown are the scaling functions O(z) (upper curve) and £(z) 
(lower) versus 2 = nt~ 1/2 at t = 16 (o), 256 (□), and 4096 
(A). For large z, both scaling functions approach the limiting 
value (4tt)- 1/2 w 0.2821. 

We can provide a relatively tight rigorous upper bound 
for the exponent ot\. We will also present a heuristic ar- 
gument, based on an uncontrolled approximation, which 
turns out to give a relatively stringent lower bound for 
a\. For both situations, our approach is based on first 
finding the ultimate survival probability of the first par- 
ticle in a finite-particle system on an infinite lattice and 
then using scaling to infer time dependence, from which 
bounds on ol\ can be inferred. Let the particles be ini- 
tially distributed on N adjacent lattice sites, with N odd. 
Ultimately, a unique particle survives which could be the 
first, the third, the fifth, etc., in the initial sequence. Let 
Si(N) be the probability that the first particle is this 
unique survivor. For N — > 00, we shall show that this 
probability scales as 

On the other hand, for t < N 2 , the finiteness of the sys- 
tem is immaterial and the survival probability of the first 
particle should coincide with S±(t) in the semi- infinite 
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system. When t becomes of the order of N, the number 
of particles remaining will be of order unity and Si(t) 
should "stick" at the value Si(N). Thus substituting t 1 / 2 
for N in Eq. (||) and equating to S\(t) gives (3\ — a%. 

First consider an upper bound for a%. A naive ap- 
proximation is to suppose that every collision between 
nearest-neighbors has the same probability to occur. For 
the initial N particles there are N — 1 collision possible 
and the probability that the first particle survives after 
the first collision is (N - 2)/(N - 1). This leaves N-2 
particles and N— 3 possible collisions, and the probability 
that the first particle survives after this second collision 
is (N — 4)/ (A — 3). The ultimate survival probability of 
the first particle in this "democratic" approximation is 



(N- 2)(AT-4)...3 x 1 
(A-l)(A-3)...4x 2 



2 



N 



(4) 



However Si(N) > 5f (N), because the first particle has 
the possibility to "escape" on its empty side, and there- 
fore collisions involving this particle are relatively less 
likely. Thus we conclude pl| 



Oil < 1/4. 



(•5) 



While we are unable to obtain a rigorous lower bound 
for oil, we have a heuristic approach that gives Si(N) ~ 
, with 7 = 3\/3/47r. This approach is based on first 
recasting the annihilation problem into an equivalent ag- 
gregation process [ fiof . In aggregation, point-like fc-mers 
perform random walks with a mass- independent diffusion 
coefficient. When two polymers of masses i and j happen 
to occupy the same site, they irreversibly aggregate into 
a heavier but still point- like polymer of mass i + j, as 
represented by the reaction scheme 
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To make the connection with annihilation, we catego- 
rize polymers according to whether their mass is odd, 
A = {A!, A s , . . .}, or even, A e = {A 2 ,A 4 , . . .}, respec- 
tively. These two classes of polymers react according to 



(7) 



In particular, the parity of odd-mass polymers is not 
influenced by even-mass polymers. Thus by consider- 
ing only odd-mass polymers, aggregation is completely 
equivalent to ARW @. 

If one associates the initial particles in ARW with 
monomers in aggregation, then the first particle survives 
in the ARW process if the mass of the leftmost poly- 
mer in the corresponding aggregation process remains 
odd throughout the evolution. With this equivalence, 



we now postulate, in the spirit of a Kirkwood approxi- 
mation, that Si(N) obeys the recursion relation 



N+2 



S 1 (N + 2)nS 1 (N)F(N + 2)& f[ F(k). (8) 



fc=i 

k odd 



Here F(k) is the probability that the first collision be- 
tween three random walks, which are initially at xo = 1, 
j/o = k — 1, and zq — k, occurs between the particles at 
k— 1 and k. While this approximation is uncontrolled, it 
turns out to give a relatively stringent upper bound for 
the true behavior of Si(N) (see Fig. (^)). 

To compute F(k), we map the problem of three ran- 
dom walks, initially at xq, yo, and zq in one dimension, 
onto a single random walk in three dimensions |l^ , pcfl . A 
collision between particles 1 and 2, and between 2 and 
3 imposes the boundary conditions that the probability 
distribution vanishes when x = y and y = z, respectively. 
This implies that the three-dimensional walk is confined 
to the wedge-shaped region defined by x < y and y < z. 
Using image techniques, the probability distribution in 
the continuum limit can be written down, from which the 
desired eventual collision probability follows after some 
tedious calculation. 

For simplicity, we give an alternative derivation which 
exploits the isomorphism between the eventual colli- 
sion probability and electrostatics p^ ]. The three- 
dimensional region x < y and y < z can be projected into 
two dimensions, with the allowed region now a wedge of 
opening angle SI = tt/3 (l^^^]. Any initial state of three 
random walks maps to a point in this two-dimensional do- 
main. To determine the co-ordinates of this initial point, 
we need to define a two-dimensional co-ordinate system 
which is perpendicular to the axis e\ = (1,1, l)/v3 
generated by the intersection of the planes x = y and 
y = z. A convenient basis is i 2 = (0, — 1, V)j\f2 and 
£3 = (—2, 1, The initial condition ro = (a;o, Uo, zq) 

has components d% = fo ■ e% = (zo — yo)/V2 and ds — 
ro ■ &3 — (—2xo + Ho + zo)/V& in the e 2 and ez basis. 
Within the two-dimensional wedge, with the horizontal 
axis defined as the locus where y = z, an arbitrary ini- 
tial condition corresponds to a horizontal displacement 
of (zo + yo — 2xo)/v / 6 and a vertical displacement of 
(zo — yo)/V2. Thus the initial point is inclined at an 
angle 



, -1 d 2 , -1 

tan — = tan 

d 3 



V3 



zo - yo 



zo + yo- 2x 



(9) 



with respect to the horizontal. 

We want to compute the probability that the random 
walk eventually hits the horizontal axis, corresponding to 
particles 2 and 3 colliding. From the isomorphism with 
electrostatics, this probability equals the integral of the 
electric field over the horizontal axis, which is generated 
by the unit charge at the initial position of the random 



3 



walk E2|. To simplify computation of this integral, per- 



form the conformal transformation w 



to open the 



wedge onto the upper half plane, so that the initial point 
is now inclined at an angle 38 with respect to the posi- 
tive real axis. In real co-ordinates, for an initial point at 
(xq, yo) with a grounded plane at y = 0, the electric field 
at (x, 0) equals 



yo 



l 



IT (x — X ) 2 + 



(10) 



The integral of this field over any interval on the x-axis 
gives the probability that a random walk which starts at 
(xo, yo) eventually hits this interval. This gives ^, where 
4> is the angular size of the interval as seen from the lo- 
cation of the charge. In our case the appropriate interval 
x = (0, oo) has angular size 4> — n — 3$, so that 

F(k) = 1 - ™ 

= 1 tan 

7T 



3V3 
2nk 1 
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k 



oo, 



(11) 



where the second line is obtained by the substitution of 
the initial condition x — 1, y — k — 1, and z = k, 
as specified by the definition of F(k). Note that this re- 
sult is easily generalizable to the case where the three 
particles have distinct diffusivities. 




FIG. 3. Simulation results for the probability that the first 
particle ultimately survives from among an initial group of 
TV adjacent particles (with TV odd) on an infinite chain. The 
lower curve represents the rigorous bound given in Eq. (4), 
and the upper curve represents the heuristic bound from 
Eqs. (8) and (11). 



Using Eq. (|j), we obtain 



N 



si(N)* n 1 



k=i 

k odd 
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N -3V3/4^ ( 12) 



Thus we arrive at the following approximate expression 
for a i 



ai = ^~ 0.2067. 

8tt 



(13) 



Let us now consider the exponent 6*2 • We write the 
survival probability of the second particle as S2 (£) = 
P(t)Sn(t), where P(t) ~ <~ 1/2 is the probability that the 
particle has not been annihilated by its single left neigh- 
bor and Sfi(t) is the probability that the particle has not 
been annihilated by any particle to its right. This latter 
probability decays as t -1 / 4 , since in the bulk the survival 
probability S(t) = S L {t)S R {t) = S R (t) 2 varies as i" 1 / 2 . 
This approximation immediately leads to ct2 — 3/4 which 
can be expected to be the lower bound. Unfortunately, 
we have been unable to construct a non-trivial upper 
bound for 0.2- A trivial upper bound, however, is pro- 
vided by the survival probability of the central particle 
in a 3-particle system. Consequently, we have the bounds 



3 3 
- ± < a 2 < -. 



(14) 



We now turn to a related and useful microscopic char- 
acterization of the reaction process, namely, the prob- 
ability that a particle is eventually annihilated by its 
n th nearest-neighbor, P(n) |23]. For homogeneous re- 
action processes, P(n) typically decays as a power law 



in n, P(n) 



, and the exponent -0 is related to the 



time dependence of the survival probability. For homo- 
geneous ARWs, for example, the particle survival prob- 
ability S(t) decaying as t~ a (with a = 1/2) implies that 
the probability that a particle is annihilated at time t is 
— dS Jf^ ~ The annihilation probabilities for given 

t and n can now be related by 



dS(t) 
dt 



dt = P{n) dn. 



(15) 



For diffusive transport, n scales as t 1 / 2 and Eq. (|l5|), to- 
gether with the defining relations for a and ip, then gives 
-0 = 1 + 2a. 

We now apply the same line of reasoning for the semi- 
infinite system. Let us define the probability that the 
first particle is annihilated by the n th particle as Pi (n) = 
n - ^ 1 ; our previous estimate for ai gives ipi = 1 + 2ai w 
1.45. Similarly for the second particle, P2M = n~^ 2 , 
with ^2 = 1 + 2a.2 ~ 2.73. As shown in Fig. H, these ex- 
pectations are consistent with our data. One additional 
interesting feature is that the second particle annihilates 
with the first particle with probability w 0.5704, while 
it annihilates with any other particle with probability 
re 0.4296. 
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the error functions. While the computation of the den- 
sity profile requires the knowledge of the two-point cor- 
relation function of the equivalent kinetic Ising model, 
the spatial distribution of the leftmost particle would re- 
quire the knowledge of all the n-point correlation func- 
tions. Although it is in principle possible to obtain such 
functions 0, this is a considerable analytical task and 
we merely use simulations to provide numerical data for 
the spatial distribution of the leftmost particle, P\ e {t(x,t) 
(Fig. ||) . This distribution obeys the expected scaling be- 
havior and is asymmetric in character, with the negative- 
z tail decaying as e _z while the positive- z tail decays as 
e~ z with z — xt~ x / 2 . From this data, we find, for ex- 
ample, that the average position of the leftmost particle 
varies as a^ieft (*) ~ 1.7 t 1 / 2 . 



FIG. 4. Simulation results for the annihilation probabili- 
ties Pi{n) (o) and P2(n) (A). For reference, the straight lines 
have slopes —1.45 and —2.73. 

Finally, we study the spatial distribution of the left- 
most particle. This refers to the extreme particle that 
currently exists and is not necessarily the initial par- 
ticle that was leftmost. To provide some perspective 
on the behavior one might anticipate, first consider this 
question for two simpler systems with the same semi- 
infinite concentration profile, c(x, t = 0) = 0(x). For 
freely diffusing particles, the long time spatial distri- 
bution in the continuum limit is the error function, 
c(x,t) = |[1 + erf(x/V4t)]. From this, the concentra- 
tion at x = remains fixed at the value 1/2, while the 
typical position of the leftmost particle is X- « —\Jt hot. 
Thus freely diffusing particles substantially penetrate the 
negative half-line. For coalescing random walks, where 
particles react by A + A — > A, the leftmost particle sim- 
ply undergoes free diffusion. Thus the concentration at 
x = vanishes as i -1 ' 2 , while the typical position of the 
leftmost particle remains at the origin. 

For ARWs, the concentration profile can be readily 
computed for arbitrary initial conditions from the direct 
correspondence to the Glauber solution of the kinetic 
Ising model || . This calculation gives 



c(n, t) 



-21 



E 

.771 = 



(-l) n I n+m (t) Jn-™(t) 



rn— 



(16) 



This profile exhibits an overall i -1 ' 2 decay of the density 
and a non-trivial spatial dependence. In the long-time 
limit, the above expression can be reduced to the scaling 
form 



;(n,t) = t~ 1/2 C(z), 



(17) 



with z = nt x / 2 and where the scaling function C(z) can 
be expressed in terms of single and double integrals of 




FIG. 5. Simulation data for the probability distribution of 
the location of the leftmost particle. Shown is the scaled 
distribution t 1 ^ 2 P\ c i t {x, t) versus scaled co-ordinate x/t 1 ^ 2 for 
t = 16 (o), 256 (□), and 4096 (A). 

In summary, we studied basic properties of a semi- 
infinite population of annihilating random walks near a 
free interface. Since particles near the interface have 
fewer potential reaction partners than bulk particles, 
these interface particles should be more long-lived than 
those in the bulk. This naive expectation turns out to 
be only partially correct. For the n th particle from the 
interface (with n = 1 corresponding to the particle at the 
interface), the survival probability S n (t) decays as t~ an , 
with a n ~ 0.225 for all odd values of n, but a n « 0.865 
for all even values of n. These exponents can be viewed 
as characterizing the surface critical behavior of ARWs 
in one dimension. 

This alternating behavior stems from the fact that an 
odd particle can eventually become the leftmost particle 
in the system and hence be long-lived. Conversely, an 
even particle will always have potential reaction partners 
on both sides and therefore has a relatively shorter life- 
time. For the odd particles, the bounds 3-^/3/8^ < a\ < 
1/4 were derived, with the lower bound non-rigorous but 
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numerically accurate, by considering the eventual sur- 
vival probability in a finite group of N particles, with N 
odd. 

The relative longevity of the interface particle is also 
reflected in the fact that the mean position of its reaction 
partner drifts slowly to the right as 1.7 t 1 / 2 . The func- 
tional form of the probability distribution of xi e ft could be 
obtained, in principle, from the n-point correlation func- 
tions of the equivalent kinetic Ising model; this appears 
to be a formidable and unenlightening task. The numer- 
ical data for F\ e { t (x,t) clearly exhibits scaling and shows 
that the position of the leftmost particle is described by 
a single length scale which varies as \ft. 
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